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Abstract. Wc consider warped five-dimensional thick braneworlds with four- 
dimensional Poincare invariancc originated from bulk scalar matter non-minimally 
coupled to gravity plus a Gauss-Bonnet term. The background field equations as 
well as the perturbed equations are investigated. A relationship between 4D and 5D 
Planck masses is studied in general terms. By imposing finitcncss of the 4D Planck 
mass and regularity of the geometry, the localization properties of the tensor modes 
of the perturbed geometry are analyzed to first order, for a wide class of solutions. In 
order to explore the gravity localization properties for this model, the normalizability 
condition for the lowest level of the tensor fluctuations is analyzed. It is found that for 
the examined class of solutions, gravity in 4 dimensions is recovered if the curvature 
invariants are regular and the 4D Planck mass is finite. It turns out that both the 
addition of the Gauss-Bonnet term and the non-minimal coupling between the scalar 
field and gravity reduce the value of the 4D Planck mass compared to its value when the 
scalar field and gravity are minimally coupled and the Gauss-Bonnet term is absent. 
The above discussed analysis depends on the explicit form of the scalar field (through 
its non- minimal coupling to gravity), making necessary the construction of explicit 
solutions in order to get results in closed form, and is illustrated with some examples 
which constitute smooth generalizations of the so-called Randall-Sundrum braneworld 
model. These solutions were obtained by making use of a detailed singular perturbation 
theory procedure with respect to the non-minimal coupling parameter between the 
scalar field and gravity, a difficult task that we managed to perform in such a way 
that all the physically meaningful conditions for the localization of gravity are fully 
satisfied. From the obtained explicit solutions we found an interesting effect: when 
we consider a non-minimally coupled scalar-tensor theory, there arise solutions for 
which the symmetries of the background geometry are not preserved by the scalar 
matter energy density distribution. In particular, the value of the "5D cosmological 
constant" of the asymptotically AdS^ space-time (which is even with respect to the 
extra coordinate) gets different contributions at -co and +oo from the asymptotic 
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values of the self-interaction potential of the scalar field. Thus, an asymmetric energy 
density distribution of scalar matter gives rise to a completely even with respect to the 
fifth coordinate space-time, in contrast to braneworld models derived from minimally 
coupled scalar-tensor theories, where both entities possess the same symmetry. 

PACS numbers: 11.25.Mj, 04.40.Nr 
1. Introduction 

Since little more than one decade an interesting alternative to the standard Kaluza- 
Klein compactification paradigm has been put to test. The kind of scenarios based 
on this alternative are known as braneworld models [I]-[6] and allow for infinite 
extra dimensions, in contrast to the Kaluza-Klein idea, where extra dimensions are 
compactified to a very small size. The mentioned alternative requires the standard 
model (SM) fields to be trapped on a 4D hypersurface, called a 3-brane. Unlike ordinary 
matter, gravitons and exotic matter are allowed to propagate in the bulk of the higher- 
dimensional manifold. Since gravity can propagate through all dimensions, the first 
important question concerning to braneworld models is to check whether they give back 
standard 4D gravity on the brane. 

In the thin braneworld models version [3]-[5] the branes are modeled by 4D delta 
functions, a mathematical complication that requires the fulfillment of the so-called 
Israel-Lanczos junction conditions, which, basically, dictate the way the brane must 
be embedded into the higher- dimensional bulk to accommodate the matter degrees of 
freedom (SM particles) that are trapped on it. These junction conditions become much 
more complicated when adding to the setup higher curvature terms, for instance [7]. 
Moreover, in these models the curvature is singular at the location of the branes, a 
drawback from the gravitational point of view that can be healed in several ways. 

On the other hand, thin brane models are just an idealization of the physical reality. 
Braneworlds, if they are to be considered as models for our world, have to be of finite 
thickness. Actually, at high enough energies, the SM particles might acquire a small 
(but non-negligible) momentum in the extra space. Indeed, the original braneworld 
idea put forth in reference [2] is consistent with a non- vanishing brane thickness ~ ''^^vk 
{mEw is the electro- weak energy scale) |j] These more realistic alternatives to thin brane 
configurations are known, generically, as thick braneworlds, or, also, domain walls. 

Thick branes might be generated in a variety of ways. One example is by replacing 
the delta functions in the action either by other distribution functions (see, for instance 
references [8] and [9]), or, alternatively, by self-interacting scalar fields minimally coupled 
to gravity [in]^[l2]. Within the framework of thin (Randall-Sundrum) braneworlds, 
bulk scalar fields have been investigated, for instance, in [131 El IIS]) where the 
Einstein-Hilbert action for Randall-Sundrum branes has been modified by considering 

I Electroweak interactions have been probed, precisely, at distances m~^^. 
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self-interacting scalar fields non-minimally coupled to gravity. In [T3] the resulting 
action has been further modified by adding a Gauss-Bonnet term, whereas cosmological 
applications have been considered in [15] . Possible influence of higher curvature terms 
in scalar-field-generated brane models has been studied in references [l6]-[26]. The 
thick braneworld approaches avoid solving the Israel-Lanczos junction conditions always 
present when studying braneworlds modeled by 4D delta functions (with or without the 
presence of scalar fields). Instead, a relevant differential equation must be solved either 
for the smooth distribution function or the scalar field of the setup. This task is, however, 
easier than facing the mathematical difficulties of solving the junction conditions when 
considering higher curvature terms in the setup. 

In dimensions larger than four (in 5D, for instance) the usual Einstein-Hilbert action 
may be supplemented with higher order curvature corrections. For some special cases, 
these corrections lead to equations of motion with at most second order derivatives 
of the metric with respect to space-time coordinates [27]. A particular combination 
containing higher curvature terms which yields second order differential equations is the 
well-known Gauss-Bonnet (GB) invariant: 

''^GB — ^ ^ABCD — 4-n: Kab + rC , [i) 

where A,B,C,D = 0,1,2,3,5. While in four dimensions the GB invariant is a 
topological term - in other words, it can be arranged into a four-divergence - which does 
not contribute to the classical equations of motion, in more than 4D the GB combination 
leads to a theory free of spin 2 ghosts due to higher derivatives and it appears in 
different higher dimensional contexts. For instance, in string theory, the heterotic string 
a correction is fixed to be the Gauss-Bonnet term in order to avoid a ghost to that 
order in the (tree-level) effective action [28]-[29]|§| However, in a second order derivative 
theory there is a sector that might contain spin 2 ghosts that are not originated by higher 
derivatives. If we wish to study braneworlds free of such ghosts, we need to carefully 
monitor the overall sign of the norm of the graviton spectrum. Within the context of 
thick braneworld scenarios generated by a self-interacting minimally coupled scalar field, 
the GB invariant has been studied in connection with the localization properties of the 
various modes of the geometry in Ref. [16] (see also, for instance, [26]). Non-minimal 
interaction of the bulk (self-interacting) scalar field has been considered, for instance, 
in Refs. [30], where the perturbative stability of the configurations was explored, and 
also within the framework of Weyl geometry in [3T] |34j. 

It seems quite natural to further explore the influence of non-minimal coupling of 
a self-interacting scalar field with gravity, within the framework of thick braneworlds 
with higher curvature terms (a Gauss-Bonnet invariant, for instance). In this work, we 
aim, precisely, at studying a 5D thick braneworld modeled by a smooth scalar domain 

§ One finds a Riemann square term in the first string tension correction and the remaining terms are 
added by hand by field redefinitions that do not afi'ect the efi'ective action to that order. However, it 
is worth noticing that if one does not choose the GB combination of the squared curvature terms, the 
corresponding graviton propagator would have terms with an arbitrary number of derivatives [29] . 
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wall non-minimally coupled to gravity with a Gauss-Bonnet term on the bulk. To be 
precise we shall consider specifically conformally fiat geometries. We shall focus here in 
the investigation of the relationship among: localization properties of the tensor zero 
mode, finiteness of 4D Planck mass and smoothness of the geometry, for a wide class 
of solutions. This will allow us to extend and generalize, in particular, previous results 
obtained within the framework of braneworlds with 4D Poincare symmetry generated 
by scalar fields minimally coupled to gravity [35] (see also [36] for a similar approach 
within braneworlds with de Sitter 3-branes). In order to do that, we analyze the tensor 
perturbations of the geometry. 

In general terms, the study of metric fiuctuations in braneworld models is 
comphcated because of the coupling between geometry and matter at the field equations 
level. On the one hand, the perturbations problem is more tractable if the background 
metric has some isometries. In this case, the metric perturbations can be classified 
according to the symmetry of the problem. If the metric respects 4D Poincare symmetry, 
then, the different perturbation modes of the geometry can be classified in a gauge- 
invariant manner into scalar, vector an tensor modes under 4D Poincare transformations. 
Furthermore, in [371 EO] it was shown in detail that each fluctuation mode evolves 
independently, then, if one studies the evolution of a speciflc mode it is not necessary to 
worry about a possible coupling with the remaining perturbation modes. On the other 
hand, in order to recover the standard 4D gravity predicted by General Relativity, the 
existence of a 4D massless spin 2 fleld localized on the brane is necessary. If there 
is a 4D massless spin 2 fleld in our model it should be described by the propagating 
tensor modes of our 5D perturbation problem. Not only the gravitational perturbation 
equations, but also the background equations are complicated and it is very difficult 
to analytically solve them. In spite of this difficulty, here we solve these equations for 
a non-singular geometry. We consider the effects of the non-minimal coupling of the 
bulk scalar fleld with the curvature as a small perturbation characterized by a small 
dimensionless parameter. Usually, in this kind of problems there are regions - called 
boundary layer regions [3H]) - where the formal expansion with respect to the small 
parameter is not valid. Then, the expansion must be redeflned for these regions. After 
having all the asymptotic expansions on the different regions, it is necessary to combine 
them to get an expansion valid on the entire domain. 

The paper has been organized as follows. In Sec. II we present the model setup and 
the corresponding fleld equations - for the respective metric ansatz - are given. In Sec. 
Ill it is further shown that 4D gravity can be localized on this particular braneworld 
and the normalization condition for a massless spin 2 fluctuation mode is derived. The 
relationship between the 4D and 5D Planck masses is obtained in Sec. IV, by integrating 
the action of the model with respect to the flfth dimension. In Sec. V a class of solutions 
for which the geometry is singularity- free is presented. Then, in Sec. VI, we further 
construct some approximate analytical solutions for the scalar fleld after flxing the form 
of the warp factor, as well as the functional dependence of the non-minimal coupling of 
the scalar fleld to gravity. Three different cases are considered by appropriate choices 
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of the free parameters. In the particular case when the Gauss-Bonnet term is switched 
off, an exact solution is obtained. Finally, in Sec. VII we present a brief summary of 
our results, and conclusions are given. 

2. The model 

Here we shall explore a thick braneworld described by the following 5D action (compare, 
for instance, with the actions investigated in references [16] and [30] ) 

s = l d'^VW\{-^ - «7^^B + ^(v^)^ - viv)}, (2) 

where a > (we choose the sign of the a coupling to be positive, although it could have 
any sign). On the other hand k, ~ where M is the 5D Planck mass, the function 

L{ip) is the coupling between the scalar field and the curvature (gravity), while V{ip) 
is the scalar field's self-interaction potential, and T^qb the 5D Gauss-Bonnet term 
defined in Eq. ([1]). The Einstein's field equations that are derivable from the action ([2]) 
are the following: 

L Rab = i^tab + VbL + ^ qab UL - e Qab, (3) 

where e = 2aK, and □ = g^^Vc^ d- The reduced energy-momentum tensor tab 
corresponding to the scalar matter content on the bulk takes the form 

2 

Tab = dA^pdsip - - Qab V{ip). (4) 

The term Qab is called Lanczos tensor and represents the corrections of the Gauss- 
Bonnet term to the Einstein equations, which can be written in the form 

Qab = 2 S'ab'T^gb ~ 2RRab + ^Rac b + ^R^ ^ Rac b d — '^RacdeRb 

Finally, the remaining terms in the right hand side of (|3]) come from the non-minimal 
coupling of the scalar field to gravity. 

The Klein-Gordon equation determining the dynamics of the scalar field is obtained 
by varying the action ([2]) with respect to ip: 

□v? + + ^ = 0, where L^ = ^. (5) 
2k OLp dip 

As in [37] let us consider a warped metric in conformally flat coordinates 

ds"^ = a'^{w) \ri^udx^dx'' — dty^] , (6) 

where the variable w is the extra coordinate and rj^y is the 4D Minkowski metric. Here, 
for simplicity, we focus in the case where the scalar field depends only on the extra- 
coordinate w. 

In terms of the above metric ansatz the field equations (l3j) and ([5]) read 

V + + ^ U" L, + L,,) -^J^n' in' + n') - [v! + ?,u') l 
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+ SW - - ^ (Sn' + 2H') = 0, 

dip K 

respectively. Here the tilde denotes derivative with respect to the extra-coordinate 
(' = while 7^ = ^, and g = L - 

The above are not independent equations. As one can easily check, if one combines 
the first two equations, the last one is obtained. 

3. Localization of gravity 

As it has been mentioned above, in linear perturbation theory there is no coupling 
between the different fluctuation modes, and hence, they evolve independently. In order 
to investigate if there is a 4D massless spin 2 field localized on the brane, it suffices to 
consider only the linear tensor fluctuations of the background metric (jS]) 



ds^ = \a?{w)-r]AB + -f^As] dx'^dx^ , where Hab = o?{w) 



2h^, 




Thus, the tensor /i^;^ is a divergence- less and trace-less rank-two tensor with respect 
to 4D Poincare transformations [8|,[37]. After taking this fact into account, perturbing 
equations ([3]) and ([5]), and performing some tedious algebra, the tensor fluctuation 
equation reads 

(q-l)' - 

2H . 

where is the 4D (Minkowski) D'Alambertian. If we redefine the field h^y as 
^ = \/ s{w)h^i, the equation ([9]) transforms into 

%y - - -D"^,. = 0, (10) 

^J S S 

where s{w) = a^q, and r{w) = a^(^q + ^^2n ) • ■'■^ order to explore the mass spectrum 
of the metric fluctuations we will assume separation of variables: ^ ^i, = ip{w)xiiu{x), 
where the field Xtiy{^) describes the 4D massive tensorial modes. Thus, equation (ITOi) 
splits into two equations: 

+ =0, (11) 

V'"-^^^V + "^'-^ = 0. (12) 

a/s S 

The function ip{w) is the w— dependent amplitude of the dynamical field Xfj.u{^) and 
defines the localization properties of the 5D field The tensor zero mode of Xiiv{x) 
can be identified with the 4D massless spin 2 field. This mode is localized on the brane 
if the associated zero mode fluctuation wave function ipo{w) is normalizable, in other 
words, the norm for ipo{w) has to be finite 

poo ^ 

(^ol^Ao) = / -ipodw <oo. (13) 
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The massless eigenstate of f ll2p is ipQ = y/s, then the above normahzation condition 
transforms into: 

/oo poo ^12 

a^qdw-Ae[a]°^^ + 8e —dw. (14) 
oo ^-oo ^ 

It is interesting to ask: what is the relationship between the 4D Planck mass and the 
normalizability condition ( !T3|) ? In the next sections we will answer this question. 

4. Planck masses 

Consider a space-time of the following form: 

ds^ = a^{w)[g^y{x)dx^dx'' — dw'^], (15) 

where gp.v{x) is an arbitrary 4D metric. Then, the relationship between 4D and 5D 
Planck masses can be obtained if we perform a dimensional reduction by integrating ([2]) 
with respect to the w coordinate. Since the coupling function L^tp) depends only on 
the extra coordinate the 5D theory ([2]) is reduced to a 4D Einstein- Hilbert effective 
action plus the corrections that come from the scalar matter and higher curvature terms 
of the bulk 

Si ~ M|i / d^x^/\^\R^ + ■ ■ • , (16) 



where the subscript 4 labels quantities computed with respect to 4D metric g^y{x). 

A way to find the Einstein-Hilbert part of the 4D effective action consists in 
considering the space-time f|T5l) as a conformal transformation of the metric qab as 
follows [391: 



Qab 9ab = a {w)gAB = a {w 



gf,u{x) 
-1 



and, rewriting the action ([2]) in terms of the quantities defined with respect to Qab- 
Since we only need to find the Einstein-Hilbert part of the 4D effective action, it is only 
necessary to apply the above conformal transformation to the terms and aT^GB 

in the 5D action. Therefore, the following expressions display the terms we need 

R = (r - 8™ - 12(W)2) , (17) 



-SRabR^^ + ^R^ - l2R{V§f - 2ARn§ + Am^^V a^ b^ 



- ASR^^Va^^b^ + 72(p^Y + 72(W)^ - 72(VaVb^9)(V^V^?9) 
+ 144(VAVBi9)(V^^9)(V^^9) + 144(W)2n^9j , 

where = \na. The quantities with a tilde are defined with respect to metric qab and 
(j2 _ C^^'^^Cabcd, where Cabcd is the Weyl tensor computed with the metric qab- 
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After separating the 4D and 5D contributions of the above quantities, substituting 
them in (|2]) and integrating with respect to the w coordinate, we obtain the following 
expression for Mpi 

~ / a%w) L{^) + —iU^ + dw (18) 

/oo 
a'^{w)qdw + ^M^e[a'\^^. 
■oo 

As one can see, the Planck mass Mpi is closely related to (t/'oIV^o); but an extra term 
arises in f lT^ . A clearer relation between Mpi and (?/'o|^/'o) can be obtained 
if we analyze the smoothness conditions for the curvature invariants. 

5. Smoothness of geometry 

A realistic thick braneworld model should not have singularities in the geometry. Since 
we are considering thick braneworlds, the corresponding warp factors, and hence the 
geometric invariants of the theory, must be smooth at the position of the branes. 
However, for some kind of solutions naked singularities can develop at the boundaries of 
the manifold [351 ESj and we should take care of the behaviour of the curvature invariants 
at spatial infinity. In consequence, let us assume that the curvature invariants and warp 
factor a are regular in the whole space-time and consider the class of solutions where 
asymptotically 

aiw — )■ oo) ~ — , (19) 
with 7 being a positive constant. For the metric ([6]) the curvature invariants are: 

R =1 [2u' + m""^ , 

R^^^''RABCD = ^,{m" + 6'H'). 
By calculating the asymptotic behaviour of the previous expressions at u; — )■ oo we get: 

R ^ w"^''-'\ R^'^Rab ^ R^^'^^'Rabcd ^ w"^''-'\ 

Thus, by imposing smoothness at infinity, the values of the constant 7 are restricted to 
the interval < 7 < 1. Under this restriction for 7, we obtain 

poo roo 

M^^c^M^ a^Ldw + --- and (^^olV^o) - / a^Ldw + ---, (20) 

J Woo J Woo 

where ■ ■ ■ denote finite terms and the interval {wooiOo) formally represents the range 
where the approximation f lT9|) is valid. Then, for this class of solutions where the 
geometry is regular, a finite 4D Planck mass implies localization of gravity on the brane. 

In order to have a physically consistent model it is not enough to satisfy the above 
conditions. Since, some terms in relations (IT^ . f ll9p and in the definition of L{ip) itself 
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have contributions with negative sign, then, it is additionally required to check whether 
the following conditions are satisfied: 

L(^)>0, (^|^)>0, M|i>0, (21) 

where a positive value of the coupling L{(f) is required in order to have a positive definite 
quadratic Hamiltonian for the metric tensor modes [30], while the positive norm of the 
graviton fiuctuation modes is required in order to have a theory free of spin 2 ghosts. 



6. Some particular solutions 

Let us consider the case where the coupling between the scalar field and gravity takes 
the form [13] (see also [30] ) 

L = l- (22) 

If the parameter ^ = the model reduces just to the one explored in Ref. [16] (the bulk 
scalar field is minimally coupled to the curvature). Hence, the parameter ^ switches 
between models with minimal and non-minimal couplings respectively. On the other 
hand, we consider a regular metric that interpolates between two asymptotically AdS^ 
[T6] space-times. This geometry is described by the following warp factor: 

ao 



(23) 



where ^ characterizes the width of the thick brane and the parameter ao is related to the 
radius of the asymptotic AdS^ space. For this geometry all of the quadratic curvature 
invariants are regular and asymptotically constant as it is shown in Fig. [H 



(Curvature) 

100 h 





-15 -10 



10 15 



Figure 1. Behavior of the different curvature invariants with respect to the extra 
space (flo = 6=1). The thick line represents the behaviour of i?^, the dashed hne 
corresponds to R^^Rab as function of x, while the thin line represents the behaviour 
of the Kretschmann scalar R^^'^'^ Rabcd- 

The scalar field (f and the self-interaction potential V{(p) can be determined using 
the first two equations of the system ([7]), since only two equations are independent. Let 
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us to replace the variable w by the dimensionless variable x = bw. In terms of this new 
variable the second equation in ([7]) can be rewritten as: 

1 + 2(1 + x^)^ (1 + x^)^ Oq (1 + x^)-^ 

In general when C, ^ it is difficult to solve the above equation. A way to 
(approximately) solve it is to assume the parameter ,^ as a small perturbation, i. e., to 
apply a perturbative analysis. It will be easier and more transparent to split the analysis 
into three separated cases corresponding to different field configurations (theories): 

I) A Gauss-Bonnet term plus a scalar field minimally coupled to gravity (^ = 0) 

m. 



II) A scalar field non-minimally coupled to gravity without the Gauss-Bonnet term 

(6 = 0). 

III) The general case where e 7^ and the small parameter ^ 7^ 0. 

Case I) Minimally coupled theory with Gauss-Bonnet term 

If one sets ^ = and considers the case where ao = 2-\/e6, the background solution 



IS 



Lp{x) = ±ipo + iff, (25) 



VT+x 



(26) 



where Vo = \ - and Vq = The ± sign describes two possible solutions of with 



8k€ ■ 

the constants (pi- In this case the 4D Planck mass does not depend explicitly on the 
parameters of the scalar field, moreover, the definition f[T^ tell us that the addition of 
a Gauss-Bonnet term to a model with a scalar field minimally coupled to gravity in the 
action slightly reduces the value of the 4D Planck mass: 

Mh ^ \^ (27) 
compared to its value when there is no Gauss-Bonnet term in the model 

Ml, (28) 

The above relation implies that the zero tensor mode is normalizable and all the 
conditions (12T]) are trivially satisfied. 

As it can be seen from ( I26p . the self-interaction potential of the scalar field V{ip) 
interpolates between two identical constant values. 

In the cases II and III we will investigate solutions of 0241) when the non-minimal 
coupling is small compared to the Gauss-Bonnet and Einstein-Hilbert contributions 
to the action ([21), characterized by the parameters e and /t, respectively. Thus, it is 
convenient to make the quantities in fl2^ dimensionless since one wishes to study the 
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field configurations independently of the choice of the units of measure. Let us consider 
the following redefinitions: ip = ^ and £ = 1^1. Under these new dimensionless 
variables the field equation ( l24l) takes the following form: 

2ex , ,,2, , 3 02 3 4e62 , , 

1 + 2 [1 + x'^)^ [l + x^Y Oq [l + x^y 

In these two cases we will assume that the solution can be formally expanded in powers 
oie% 

(t) = (t)Q{x)+e<j)i{x)+o{e<j)i{x)). (30) 

By substituting this expansion into ( l29l) one can find approximate solutions of the 
field equation. However, a subtlety arises when we try to apply this naive procedure 
to (l29ll since in this case the term containing the second derivative is multiplied by 
the small perturbation parameter, namely, if one sets e = in f l29|) the differential 
equation that arises is no longer of second order, but of first order. Therefore, it is 
not possible to generate all the approximate solutions to the field equation ( 129|) with 
two arbitrary boundary conditions for at x = ±oo, because the differential equations 
for 00 and 0i will be of first order. Thus, the study of the boundary value problem 
poses a dilemma: there is only one unknown arbitrary constant but there are two fixed 
boundary conditions to be satisfied. This leads to what is generally known as a singular 
perturbation or boundary layer problem [381 SOI SH 1121 US]. 

In this kind of problems, in general, it is not possible to find a single uniformly 
valid asymptotic expansion for the field, in agreement with two arbitrary boundary 
conditions. What one can do is to propose two different asymptotic expansions for 
regions that contain the boundaries x = —oc and x = oo, such that each one of them 
is uniformly valid in the corresponding region. Although, as we will see, in our case it 
is not enough to have two asymptotic expansions: we shall require another asymptotic 
expansion in the neighborhood of the origin (x = 0). 

In addition to the mentioned complications, in this kind of problems usually there 
is a region of rapid variation of the field and/or of its first derivatives, known as the 
boundary layer. In this region the term that contains the second derivative in the field 
equation is no longer negligible, thus the expansion (l30l) and the corresponding solutions 
are not valid in this region. 

In order to find an asymptotic expansion uniformly valid in a boundary layer located 
at Xb, characterized in size by the function 5{e), it is suitable to magnify this region by 
rescaling the variable x with the aid of the stretched or boundary layer variable: 

( = with 5{e) = o(l) when e ^ 0. (31) 

II If the first-order term £(/)i(x) in (|30p is uniformly "small" with respect to the zeroth-order (t)o{x) 
when £ — > over some region in the domain of variation (2?) of 0, in other words, if e(j)i{x) ~ o{(j)Q{x)) 
when £ ^ over some region € I?, the expansion is called asymptotic expansion up to first order of the 
field on this region. 
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In terms of this variable the field transforms into: 

0(x) = 0(x, + 5(£)C) = $(C). (32) 

The next step is to consider a different expansion for $(C) that we hope to be uniformly 
valid on the boundary layer region. How to choose 6{e) in order to satisfy this condition 
will be discussed later. 

On the other hand, in order to obtain an approximation on the whole domain T>, it 
is necessary to know how to match two different adjacent asymptotic expansions. The 
principal hypothesis of this method is to assume that there is an intermediate region 
where two different asymptotic expansions give the same result. By following [IQ], the 
idea consists in defining on the intermediate domain a new stretched variable: 

Co = r r \ J with 6o{e) = o(l), 6{e) = o{6o{e)) when e — )■ 0, (33) 



and rewriting the asymptotic expansions in terms of Co up to some ordeinJ and then, 
after expanding both of them, asking for their equality in terms of an arbitrary function 
6o{e) under the restrictions written in fl33l) . This matching procedure serves us to fix 
the constants that appear in the expansion for $(C) on regions which are not connected 
with the boundary conditions. 

A description of the solution on certain regions consists of two expansions which 
must be combined to form a composite expansion. This is done by adding the expansions 
and then subtracting the part that is common to both, yielding an approximation to 
the solution on the above mentioned region [13]. After determining all the composite 
expansions uniformly valid on different regions, the resulting approximation to the 
solution valid on the whole interval can be obtained by joining together all of them. 

Case II) Non-minimally coupled theory without Gauss-Bonnet term 

In this case the equation (l29l) can be written as 



'"^ ^^' + <P'\^-^)-l^7V^ = -7r^^ (34) 



1 + X2^^ ^ ' '2 (1 + X2)2 {l+X^y 

and the consistency conditions (121 p can be reduced to the following pair of constraints: 
L(0) = 1 - > and M^j > 0. 

Ila) Case without boundary conditions 

These last restrictions do not require to fix any initial or boundary conditions. 
Thus, to begin with, let us study the class of approximate solutions generated by (150]) . 
By substituting this expansion into (l34l) it is not difficult to obtain the equations for 
the first and second approximations, (poi^x) and 0i(a;), respectively: 

^? - = 0. (35) 



^ Not necessarily the same order for both asymptotic expansions. 
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The above equations give us the following solutions 

arctan(x) + 



^0 ; 



±^arctau(.) 

2 ^ ^ 



V3 



' arctan(x) =F — — arctan 



3 

— -Aq arctan2(x) + Af, 



(37) 



(38) 



where the signs ± mean two possible solutions for each order of approximation, and 
and Af are arbitrary constants. One can check that the above solutions generate an 
asymptotic expansion 0^ = 0q + e(j)f on the whole domain V = (— oo, oo) of the field 
variable. In Fig. |2]it is shown that there are configurations of the scalar field with a kink- 
like behaviour. Moreover, in contrast to the previous case, the self-interaction potential, 
in general, interpolates between two different negative constant values as shows the Fig. 
131 these asymptotic values can be written as follows: 



V+(±oo) 



V^"(±oo) 



3b^ 



e 



-8 + £ (2A^ T VStt 



where and V~ correspond to the solutions 0+ and (f)~ respectively. This asymmetric 
asymptotic behaviour of the matter energy density of the system may seem surprising 
if one takes into account the fact that the space-time is asymptotically AdS^. However, 
the uneven character of the energy distribution of the scalar matter at — oo and +oo is 
compensated by its non-minimal coupling to gravity, rendering an asymptotically AdS^ 
space-time which is even with respect to the fifth coordinate. Furthermore, up to first 
order in e the 4D Planck mass is 



-TT 



(39) 



As we expected from f|T9l) . the above formula tells us that the effect of the non- minimal 
coupling between the scalar field and gravity reduces the value of the 4D Planck mass 
compared to its value when there is only minimal coupling. By taking into account the 
requirements L{ip) > and M^j > mentioned above, one can show that the values of 
e are bounded 

2 



e < 



(40) 



In Fig. H] it is shown the function L{(f)) = 1 — for several sets of values of the 
parameters. As one can see, the condition L{(f)) > is obviously satisfied. 

As one can see from ( 137]) and fl38l) . after fixing the sign of the solutions, only one 
free boundary condition is left. In consequence, it is not possible to generate all of the 
approximate solutions of f l5^ by using only one expansion of (p. 

For arbitrary boundary conditions it is difficult to solve (!34|) . Hence, here we 
solve the field equation just for the case where the field vanishes at x = ±oo. These 
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Figure 2. Graph of the function 4){x) up to first order in e for case Ila). In the figure 
we set £ = 0.01 and = Af ~ 0. The thin fine represents 0+ = + Ef^i^ and the 
thick one describes (/)~ = c/jq" + E(\\ . 
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Figure 3. Self-interaction potential Vij^ix)) up to first order in e for case Ila) (in all 
cases we have chosen k = 1, ap = 1, 6 = 1 and e = 0.01). The following values of the 
constant parameters have been chosen: for = - thin line, = 3 - thick 

line, while the dashed line corresponds to the choice = 3, for 0^. 

boundary conditions are interesting because these are the only boundary conditions that 
will appear in the next relevant case where ^ 7^ and e 7^ 0. 

lib) Case with imposed boundary conditions 

Let us consider the equation (IMj) under the following boundary conditions: 

0(-cx)) = 0(cx)) = 0. (41) 

First of all, it is easy to show that if 0(x) is a solution of (134|) under ( 14T|) . then, 0(— x) 
is a solution of the field equation too, with the same boundary conditions. This implies 
that the solution of the field equations under the conditions ( HTl) is an even function. 

Let us solve perturbatively (IMjl under the boundary conditions (j^Tj) . The expansion 
outside the boundary layers is described by the solutions ( l37|) and ( !38l) . As one can 
see, it is not possible to construct a single expansion which satisfies both boundary 
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Figure 4. Behaviour of the function L{4)) - up to first order in e - vs the variable x 
for the case Ila). We have chosen three different sets of values of the free constants 
(e ~ 0.01 in all cases): i) Aq = 0, 0^ - thin line, ii) Aq = 3, 4'q (thick line), and iii) 
Aq = 3, 0q - dashed line. 

conditions at the same time. Therefore, we need to find asymptotic expansions vahd on 
the boundaries. In general, this is a difficuh task, but in our case we shall see that if 
one chooses any expansion or for a domain that contains the boundary x = +oo 
and the remaining solution for a region that contains the point x = — oo, the constants 
that appear in the solutions can be fixed in such a way that the boundary conditions 
can be completely satisfied. 

In what follows, we will denote with the index a the case where 0~ is defined on 
the region that contains the boundary x = —oo and (p^ on the region that contains the 
boundary x = +oo, otherwise the subscript (3 will be used. To begin with, let us discuss 
in detail the case a. The other case is similar, therefore, only the final results will be 
presented. 

By setting = = Aq = -^n and Af = A]; = Ai = - ij, we can 

write the approximate solution to flM|) under fHTj) for regions that contain the boundaries 
in the form: 

0~ = — -\/3 arctan(x) + Aq ^£:[2 arctan(x) — arctan(a;) — arctan^(a;) 

+ VSAo arctan^(x)] + eAi, for allx G V~, (42) 
0^ = v^arctan(x) + Aq + ^p£:[2 arctan(x) — y4Qarctan(x) — arctan^(x) 

- VsAo arctan^(x)] + eAi, for allx G V^, (43) 

where the domains V~ and contain the points x = — oo and x = +oo, respectively. 

If one naively supposes that there is no intermediate boundary layers between the 
two boundary regions, the next step is to find a common region where the two expansions 
are valid. In our case, it is not difficult to find it because and 0^ can be matched 
order by order at x = 0. Thus, the domains can be taken as P~ = (— oo,0] and 



Thick braneworlds generated by a non-minimally coupled scalar field and a Gauss-Bonnet termlG 



T)^ = [0,00), respectively. These two domains allow us to define a general solution 
(valid on the entire domain) for the field. As one can see from (H2l) and P3|) . the latter 
solution is continuous but its derivative is not continuous at x = 0. In other words, 
on a neighborhood of the origin the first derivative of the field changes "rapidly" from 
negative to positive values. Therefore we have a boundary layer behaviour of the first 
derivative of the field on a neighborhood of the point x = 0. 

The above analysis tell us that our initial hypothesis about the absence of 
intermediate boundary layers is false, then, we need to find an asymptotic expansion 
valid on a neighborhood of the origin. As we mentioned before, it is convenient to define 
a new variable 

X 



c 



with 6{e) = 0(1), when 5 — )■ 0, 



6{ey 

and to perform a new expansion for 
origin a; = 0: 

In terms of the new variable (l44l) the equation ([3 

c 



(44) 



$(C) over the domain that contains the 



d<^ 

{scy 'dc 



can be rewritten as: 
2 



2 I 



2^2 



(45) 



(46) 



By substituting (145!) into the previous equation, we get 



dC 

lie 



+ e 



+ 2 



$0^ 

° de 

d^o d^i 
dC dC 



$1 







dC dC 








dC ' 





dC 



e6' 



2(%- 



d^r 



dC 



(<fc 



-35' + 



+ 



(47) 



In order to determine $0 and $1 we need to know how to choose 5{e). The idea is to 
make a well-balanced choice of the variable ( or, equivalently, of the S{e) function, such 
that the equations for the approximations $0 and $1 contain as much information as 
possible as e — 0. This selection is called the distinguished limit of the expansion (H5l) 
[38]^|l2]. Let us put it in different words: if we choose 6{e) in a way different from 
the distinguished limit, all the terms that arise in the equations for $0 and $1 will be 
contained in the equations for these functions in the distinguished limit. 

In our case we can show that 6 = e corresponds to the distinguished limit of ( H5l) . 
therefore, by using ( Wf} the equations for the zero and first order approximations of $ 
on the boundary layer region are: 

d% 



dC 



d'^i fd^i 



dC 



0, 



-3. 



(48) 
(49) 
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By solving the above system of equations, the field $(C) on the boundary layer domain 
adopts the form: 



$(0 = $o + £ 



C2 — $0 In cosh I — (C + Ci) 







(50) 



where Ci, C2 and $o cire constants. These constants can be determined by matching the 
expansion of f l42p . up to some order, with the above solution on a negative neighborhood 
of the point x = 0, and by matching the expansion of the solution f H3|) . up to some order, 
with (ISO]) on a positive neighborhood of the point x = 0, respectively. Let us next define 
the intermediate stretched variable: 

X 

Co = T^T' '^i^h 5Q{e) = o(l), e = o{So{e)) when e ^ 0. (51) 

By expanding the zeroth order approximations of ( 142|) and ( l43l) and the boundary layer 
solution (!50|) in terms of 6o{e), we gelE|: 



±y3<5oCo - V3 



TT 



$± ~ $0 ± V36oCo + e(±y3ci + C2 + $oln2), (52) 

where in the ± symbol, the — sign stands for quantities defined on a negative 
neighborhood of the origin, similarly, the + sign denotes quantities defined on a positive 
neighborhood of the origin. 

By matching the above expansions, the following values for the constants $o; Ci 
and C2 are obtained 

$0 = -V3^, 
ci = 0, 

C2 = VS^ In 2. 

Notice that, up to this point, the description of our approximate solution is realized 
by three well-defined asymptotic expansions (which form two pairs), that can be joined 
together to get a solution that is valid on the entire domain V = {— oo < x < +oo}. One 
pair of these asymptotic expansions is uniformly valid on the region {x < 0} G V, while 
the other pair is uniformly valid on the region {x > 0} G V. Thus, we can construct 
two composite expansions defined on the regions x < and x > 0, respectively, that 
will be collected to generate an approximate solution defined on the whole domain of 
variation V. 

As we mentioned above, the composite expansion for the region x < (or x > 0) 
is obtained by adding two asymptotic expansions uniformly valid on the corresponding 
region and expressed in terms of the same variable, and then by subtracting the part 
that is common to both of them. In our case it is convenient to construct the composite 
expansions for the positive and negative regions (x < and x > 0) separately. 

+ In order to include the first order approximation 0^0,(2;) in the matching process, it is necessary to 
calculate the contribution of £^$2(C) to ttis expansion (j45|) . However, this ease is analytically more 
complex and, thus, will not be considered here. 
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In other words, the composite expansion (either for x < or x > 0) can be expressed 
as follows 

= + * (^) - 0co™(x), (53) 

where 4>com{x) is the common part to both asymptotic expansions. It should be noticed 
that in our case, by virtue of fl52l) . the common functions for the expansions on the 
X < and x > regions can be described in a single functional form 

(pcom{x) = VS ^|x| - . (54) 




Figure 5. The graph (a) shows the solutions for the field (j) in the case lib); the thin 
line represents 0" (x) and the thick line represents (f>^ (x) . Both profiles asymptotically 
tend to zero. In graph (b) we represent the profile of (j)^{x) in the neighborhood of the 
point a; = 0. In these figures we set e = 0.01. 



V 




V 




Figure 6. The graph (a) shows the self- interaction potential vs x for the case lib). 
In graph (b) we zoomed the self-interaction potential in the neighborhood of the point 
a; = 0. In these figures we set k = 1, ao = 1, 6 = 1 and e = 0.01. 



Thus, the solution is: 



0°(x) = 0(x) = Vs < arctan |x| — ^ 







7r\ 


X — 


X - 











TT 

2' 



2cosh 



(55) 
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Figure 7. The graph (a) shows the function £(0) vs x for the case lib). In graph 
(b) wc zoomed this function in the neighborhood of the point a; = 0. In these figures 
we set e ~ 0.01. 



By following the same scheme described above, we can show that </)"(x) = —(()'^{x). 

As it is shown in Fig. E^, both solutions are bounded, continuous even functions, 
with continuous derivatives on the whole domain V. In Fig. [5)d we have zoomed (j)^{x) 
on a region near to the origin of coordinates to appreciate more clearly the smooth 
behaviour of this profile. Furthermore, like in the case I, the self-interaction potential 
interpolates between two identical negative constant values given by (see Fig. [6^): 

V(oo) ~ 



On the other hand, the 4D Planck mass is 



M|i~ ^^{2-3£(7r-2)}, (56) 

again, the non-minimal coupling effects diminishes M|j with respect to a simpler model 
where the field is minimally coupled to gravity and the Gauss-Bonnet term is absent. 
By using L{ip) > and the above obtained result fl56|) . the restriction for the values of 
.is g 

In Fig. [Tti it is shown the function L{(j)) for both solutions. As we wished, the 
condition L{(f)) > is satisfied. 

Case III) Non-minimal coupling and Gauss-Bonnet term 

Like in the first case, the relation ao = is imposed, then, the field equation 

(129|) is transformed into: 



1 + ' ^ V . 2 (l + x2)2 (l + x2)3' 

In this case we shall modify the condition (^/;|?/') > of fl^Tl) by a stronger and more 
tractable one, where the weight function rjs is not negative. Under this redefinition of 
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(I2T]) the admissible profiles for are constrained by the following condition: 

m = IT^ - ^ 0' (58) 

then, when x — > oo the field vanishes. Based on this fact, we will study approximate 
solutions to ( l57j) under the boundary conditions: 

0(-oo) = 0(oo) = 0. (59) 

Of course, this boundary value problem is weaker than ( 157|) under the restrictions 
( I2TI) and ( |58l) . Thus, after finding the solutions, it is necessary to check whether such 
restrictions are fulfilled. 

As mentioned above, we will find solutions of ( l57|l in terms of the expansion ( 130|1 . 
By substituting such an expansion into ( 1571) . the equations for the zeroth and first 
approximations become: 

^ - 0, 



W-<Po-'Po0o-Y^ + 2(rT^ = O' 



and possess the following solutions: 



- ^^^^/T^-^o' 



± 



0± = ^ |-|S= + 3 [5 + [B^r] arcsinh(x) - 4v^5o'= ln(l + x^)} + 

where and Bf are constants. As in the previous case, for each order of perturbation 
we have two solutions characterized by the signs ±. Since we have well-defined boundary 
conditions, we must set the values of the above constants. Again, with a single solution 
it is not possible to satisfy both conditions simultaneously. Then, motivated by the 
case II we have two possibilities of finding asymptotic expansions which are valid on 
the boundaries. A first one where 0~ is valid on a region that contains x = —00, and 
0"'" is valid on a region that contains x = +00. By following the same notation used 
in the case II, we will denote this choice by an a index. The second one corresponds 
to choosing the expansions in reverse order and will be denoted by a /3 index. Let us 
consider the first case in detail. 

Like in the case II, a boundary layer behaviour of the first derivative on the 
neighborhood of the origin arises. In order to overcome this problem we repeat the 
same procedure applied to the case II: we perform the change of variable (jH]) and 
express ( l57|) in terms of it: 



(1 + (5C)¥' 

where 0(C) is the field written in terms of the ( variable. When one compares ( 147|) to the 
above equation, the only difference appears in the right hand side of the equations, but. 
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after expanding (15U]) on the boundary layer region, we only consider the first-order term 
in S. Thus, the equations for the two first approximations of the field on the boundary 
layer domain are the same that we obtained in the previous case and the solution takes 
the form: 



e(C) = 00 + £ < c/2 - ©0 In cosh 



V3 
$0^ 



(61) 



where 0o is constant and describes the zeroth order approximation of the field and di, 
d2 are arbitrary constants. By matching the three expansions we have: 

60 = -V3, 

di = 0, 

d2 = 73 In 2. 




Figure 8. The solutions for the field </> in the case III; the thin line represents 

and the thick line denotes (f)^{x). Both profiles tend asymptotically to zero. In the 



figure we chose £ = 0.01. 

Finally, the solution for the field is 

0"(x) = 4){x) = V?,[^^^=-l - |a;|+eln[2coshf-)l 1. (62) 

By applying the same procedure as in case lib), one can show that (j)^{x) = — 0"(a;). 
From Fig. |H]we see that both profiles of the field are continuous, bounded even functions, 
with continuous derivatives on the whole domain. The asymptotic value of the self- 
interaction potential reads: 

and its profile is shown in Fig. |9l it is similar to the one illustrated in Fig. [6^. 
Furthermore, the value of Mpj is 



M|,~^ (--.). (63) 



M3 /4 
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Figure 9. Self-interaction potential vs x for the case III (ac = 1, oq = 1, 6 = 1 and 
e = 0.01). 

By comparing the previous expression to f l27p one can see the effect of the non-minimal 
couphng: it reduces the value of the 4D Planck mass compared to its value in a model 
in which the scalar field is minimally coupled to gravity and there is a Gauss-Bonnet 
term on the bulk. The consistency conditions fl2T]) together with the modification fl58|) 
imply that 

e < I (64) 

As one can see from Fig. [101 the consistency condition f l58l) is fully satisfied for these 
constrained values of e. 

f(x) 




-6 -4 -2 2 4 6 



Figure 10. The figure shows the fulfillment of the consistency condition ([58| . As 
before we choose e = 0.01. 

In order to compare the cases II and III we plot (f)^{x) and the self- interaction 
potentials for both cases in Figs. [H] and [121 respectively. Although there are some 
small differences for these quantities near the origin of the extra coordinate, similarly to 
the results obtained in [H], one can see that there is no significant difference between 
both cases. 
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Figure 11. Comparison of the solutions. The thin hne represents (f)^{x) for the case 
lib) and the thick hne represents for the case III. In both profiles we have 

arbitrarily set e = 0.01. 

V 




Figure 12. Comparison of the self-interaction potentials. The thin line represents 
the self-interaction potential for the case lib) and the thick line represents the self- 
interaction potential for the case III (k = 1, oq = 1, 6 = 1 and e ~ 0.01). 



As we have already mentioned, in all the analyzed cases the consistency conditions 
(pTll are satisfied. Since the couplings L{(f)) are bounded, the corresponding zero 
tensor modes are localized on the brane. Therefore, in our model ([2]) there are field 
configurations that give rise to a regular and asymptotically AdS^ geometry with a 4D 
massless spin 2 field localized on the brane. Unlike the zero tensor modes, it is much 
more difficult to analytically find the full massive spectrum. In [30] it was shown that in 
the case II, where the Gauss-Bonnet term is absent, the massive spectrum is continuous 
for the class of solutions (fT9l) . In other words, for the case II we have a localized zero 
tensor mode and a tower of continuous massive modes without a mass gap. It is more 
difficult to characterize the massive spectrum for cases I and II because the equation for 
the mass spectrum does not have the form of a Schrodinger equation, at least in terms 
of the coordinate x = bw; thus, these cases will not be considered here, but they will be 
treated in a later work. 
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7. Conclusions 

We have explored a thick braneworld model where the matter scalar field is coupled 
non-minimally to the Einstein-Hilbert term, in addition to this, there is a Gauss-Bonnet 
term in the bulk. We compute the 4D Planck mass in terms of the 5D Planck mass and 
quantities related both to the matter and the geometry. In contrast to theories where 
the matter field is minimally coupled to the Einstein-Hilbert term, in our model Mpi 
depends explicitly on the matter content of the bulk due to the non-minimal coupling 
of the scalar field to gravity. In addition to this, by imposing certain natural conditions 
on the parameters of the model (see ( 12T1) ). the expression for the 4D Planck mass (see 
(fT9|) ) tells us that if a non-minimally coupled scalar field and/or a bulk Gauss-Bonnet 
term are/is considered, the predicted 4D Planck mass will be smaller than that resulting 
in a model where the scalar field is minimally coupled to gravity and the Gauss-Bonnet 
term is absent. 

On the other hand, a relation among smoothness of geometry, finiteness of the 4D 
Planck mass and localization of the tensorial modes was studied for a wide class of 
solutions. Our results show that if the geometry is regular, a finite 4D Planck mass 
implies localization of gravity on the brane. In the general analysis described above 
some assumptions were made about the regularity of the geometry and its asymptotic 
behaviour. 

We further explored an example in which we applied this general analysis to a 
regular, asymptotically AdS^ geometry. Since the 4D Planck mass and the normalization 
condition explicitly depend on the profile of the scalar field, it is important to solve 
the equation for and then to check carefully whether the conditions (|2T|) are fully 
satisfied, since in such quantities there are terms with negative sign contributions. 
We perturbatively solved the equation for ip by considering a small parameter defined 
in terms of the strength of the non-minimal coupling. In order to solve the latter 
equation it was necessary to apply the singular perturbation method since the scalar 
field appeared to have different scales of variation on different regions, a situation known 
in the literature as the boundary layer problem. The application of such a method is 
a difficult task in general, since the involved differential equations are hard to solve 
analytically in exact or approximate form. However, for a wide class of solutions we 
managed to get perturbative analytic expressions for ip which satisfy all the physically 
meaningful conditions ( 12T]) . 

In order to elucidate the physical effects of the inclusion of the non-minimal coupling 
and of the Gauss-Bonnet term, we further considered three cases in which we switched 
on/off the corresponding parameters. When the model has minimal coupling with a 
Gauss-Bonnet term, the solution can be obtained exactly and ip has the form of a 
kink/ anti-kink. In the second case, when the non-minimal coupling is considered and the 
Gauss-Bonnet term is switched off, there was obtained a class of perturbative solutions 
with kink-like behaviour. 

For this latter case we can observe an interesting effect: there are some 
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configurations where the self-interaction potential of the scalar field approaches different 
asymptotic values at — oo and +00, mimicking two distinct cosmological constants at 
both ends of the extra coordinate. Since the space-time is asymptotically AdS^, this 
means that the non-minimal coupling of the scalar field to gravity also contributes to 
the "total" cosmological constants at —00 and -|-oo, compensating the unevenness of 
the scalar energy distribution in order to asymptotically render an even AdS^ space- 
time. Usually, in General Relativity the symmetries of the geometric background are 
preserved by the matter energy distribution as a consequence of self-consistency of the 
Einstein equations. We see that this is not the case anymore when one considers a 
system with non-minimal coupling between the scalar field and gravity. In particular, 
for a symmetric or even with respect to the extra coordinate geometry, the "total" 
cosmological constants of the asymptotically AdS^ space-time get distinct contributions 
from the asymptotic values of the self- interaction potential of the scalar field at —00 
and -|-oo, and the non- minimal coupling between the scalar matter and gravity. 

We further imposed asymptotically vanishing boundary conditions on the field 
configuration of the second case, yielding a physically meaningful solution which is 
even with respect to the fifth coordinate. This solution has been compared with the 
solution of the third case in which the Gauss-Bonnet term is also considered. We found 
that both solutions have a very similar behaviour. 

Moreover, for these two last cases we also found that the consistency conditions 
fl2T|) impose bounds on the perturbation parameter e and that they render a smaller 
4D Planck mass compared to the case in which both the non-minimal coupling and the 
Gauss-Bonnet term are absent. 

It will be interesting to further consider non-minimal coupling of the scalar field 
with the Gauss-Bonnet term. The investigation of the resulting model will be the subject 
of forthcoming work. 
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